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ABSTRACT 

We combine new ASTE/AzTEC 1.1 mm maps of the galaxy pair NGC 
1512/1510 with archival Spitzer IRAC and MIPS images covering the wavelength 
range 3.6-160 fim from the SINGS project. The availability of the 1.1 mm map 
enables us to measure the long-wavelength tail of the dust emission in each 
galaxy, and in sub-galactic regions in NGC 1512, and to derive accurate dust 
masses. The two galaxies form a pair consisting of a large, high-metallicity spi- 
ral (NGC 1512) and a low metallicity, blue compact dwarf (NGC 1510), which 
we use to compare similarities and contrast differences. Using the models of 



Draine fc Lil (120071 ). the derived total dust masses are (2.4±0.6) x 10^ Mq and 
(1.7±3.6)xl0^ M0 for NGC 1512 and NGC 1510, respectively. The derived ratio 
of dust mass to H I gas mass for the galaxy pair, Md/Mni ~0.0034, is much lower 



Draine et al. 



(by at least a factor of 3) than expected, as previously found by 
(120071 ). In contrast, regions within NGC 1512, specifically the central region and 
the arms, do not show such unusually low Md/Mni ratios; furthermore, the dust- 
to-gas ratio is within expectations for NGC 1510. These results suggest that 
a fraction of the H I included in the determination of the Md/Mni ratio of the 
NGC 1512/NGC 1510 pair is not associated with the star forming disks /regions 
of eith er galaxy. Using the dust masses derived from the models of iDraine &: Li 
(120071 ) as reference, we perform simple two-temperature modified-blackbody fits 
to the far-infrared/mm data of the two galaxies and the sub-regions of NGC 
1512, in order to derive and compare the dust masses associated with warm and 
cool dust temperature components. As generally expected, the warm dust tem- 
perature of the low-metaUicity, low-mass NGC 1510 (T„~36 K) is substantially 
higher than the corresponding warm temperature of the high-metallicity spiral 
NGC 1512 (T„~24 K). In both galaxies (albeit with a large uncertainty for NGC 
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1510), our fits indicate that a substantial fraction (>93%) of the total dust mass 
is in a cool dust component, with temperatures ~14-16 K for NGC 1512 and 
~15-24 K for NGC 1510. This result is similar to what is determined for a few 
other nearby galaxies. In contrast, the warm dust component in the sub-galactic 
regions of NGC 1512 represents a much larger fraction of the total dust content, 
in agreement with the fact that all three regions have higher specific star forma- 
tion rates than the average in the galaxy; in the center, the warm dust represents 
about 40% of the total, while in the arms the fractions are close to ~20%. 

Subject headings: galaxies: individual (NGC 1510, NGC 1512); galaxies: ISM; 
galaxies: starburst; (ISM:) dust, extinction 
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INTRODUCTION 



Interstellar dust provides a source of attenuation of the light from galaxies, but also a 
source o f information on th eir metal content, and, possibly, a tracer for the molecular gas 



content (jPraine et al 



20071 ). On average, about 50% of the UV-optical light from galaxies 



i s re-processed by their own dust into the mid/far-infrared and mm wavelength range 



(jPole et al. 



20061). 



Determining accurate dust masses for galaxies has been, traditionally, a difficult step 
since measurements at infrared wavelengths shorter than ~100-200 /im (as accomplished 
by IRAS, ISO, and the Spitzer Space Telescope) are sensitive to the inferred dust mean 
temperature(s) and the adopted dust emissivity. Wavelengths longer than a few hundred 
yum probe the long wavelength tail of the modified Planck function for typical galaxy dust 
emission, and the inferred dust masses are, therefore, less sensitive to the dust temperature. 
Long-wavelength measurements also offer additional leverage for constraining, when used 
together with shorter wavelength infrared data, the dust temperature itself. 

Recent applications of sub-mm data used to constrain the properties of the dust 
emission from galaxies include a nearby sample of about 100 galaxies observed at 450 /im 
and 850 um with SCUBA (The SCUBA Local Universe Galaxy Survey, a.k.a., SL UGS, 



Dunne et al. 



2000 



Dunne fc Eales 



2001 



Seaquist et al.il2004t iVlahakis et al. 



20051 ). and 



an addi tional ~20 galax i es fro m the SINGS sample (the Spitzer Infrared Nearby Galaxies 



Survey, 



Kennicutt et al 



20031 ). where Spitz er data in the ran 



combined with SCUBA 850 /im detections (IDraine et al 



ge 3-160 /im have been 



20071 ). More recently, sub-mm 



data have been comb ined with Spitzer data for 11 nearby galaxies from the SLUGS sample 



( IWillmer et al. 



20091 ). thus increasing the total number of galaxies for which the dust masses 



are measured accurately from multi-wavelength spectral energy distributions, spanning 
from the mid-IR to the sub-mm. The SLUGS galaxies tend to have, on average, colder 
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dust temperatures th an the galaxies 



part), as estimated by 



Willmer et al 



in the SINGS sample (of which NGC 1512/1510 is 



(|2009[ ) on the basis of the far-infrared/sub-mm colors. 



The Spitzer data on nearby galaxies offer the unique opportunity to analyze the dust 
emission characteristics not only of whole galaxies, but also within sub-galactic regions, 
thanks to its angular resolution in the range 2"-38" for the 3-160 fim range. We combine 
the SINGS data on the interacting galaxy pair NGC 1512/1510 with data at 1.1 mm 
wavelength from ASTE/AzTEC, which has 28" resolution, to investigate the infrared/mm 
spectral energy distribution (SED) of dust for those two galaxies, and for regions within the 
larger NGC 1512. Our goal is to derive accurate dust masses for the two galaxies, and to 
compare the SEDs between sub-regions. 

The interacting pair NGC 1512/1510 is of particular interest because its two member 
galaxies are remarkably distinct in nature: the primary galaxy, NGC 1512, is a large, 
metal-rich barred spiral, while its companion, NGC 1510, is a low metallicity blue compact 
dwarf (BCD) galaxy (Tabled]) which is expected to have less ability of self-shielding for the 
dust. Both galaxies are characterized by intense starburst activity in their central regions, 
and do not show evidence for presence of AGNs. This pair is thus ideal for investigations of 
dust SEDs, and their dependence on each galaxy's properties. 

The paper is organized as follows: section 2 describes the data and, for the 
ASTE/AzTEC data, provides some detail on the observations and data reduction; section 
3 is devoted to the data analysis, including the photometry, and photometric corrections 
and uncertainties; section 4 presents the dust model fitting and results; the discussion and 
conclusions are in section 5. 
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DATA 



2.1. IRAC and MIPS Data 



Spitzer maps with both the IRAC (3.6, 4.5, 5.8, and 8.0 fim) and MIPS (24, 70, 
and 160 fim) are available for both galaxies, through the high le vel data products of t he 



Kennicutt et al 



20031 ) 



SINGS Legacy project (Spitzer Infrared Nearby Galaxy Survey, 
Although the SINGS observations targeted NGC 1512 only, its companion galaxy is 
serendipitously present and detected in the same maps. The SINGS ob servation strateg y. 



data r educ tion procedures, a nd map sensitivity limits are described in 



mom and 



Dale et al 



Kennicutt et al. 



( 120051 ). The angular resolution is ~2" for the four IRAC bands and 



6", 17" and 38" for the MIPS 24, 70 and 160 /xm bands, respectively (Figure [I]). 



2.2. AzTEC/ASTE Data 



AzTEC is a 144 e lement bolometer ar ray currently configured to observe in the 1.1 mm 



atmospheric window (iWilson et al. 



2008a] ). Observations of NGC 151 2 were made with 



AzT 



]C on the Atacama Submillimeter Telescope Experii nent (ASTE, 



20081 ) using a network observation system N-C0SM0S3 (IKamazaki et al. 



Ezawa et al. 



2004, 



20051), which 



provided a 28" circular main beam. During the period 7-11 October 2007 w e performed 



Wilson et al. 



(l2008bl ) centered on 



25 lissajous observations identical to those described in 
(04^03"54^28, -43°20'55.9"). The coaddition of the 25 observations results in a roughly 
circular map with uniform coverage over the central 18' diameter with a depth of ^1.5 mJy. 
The structure of NGC 1512 is well resolved, with a maximum signal-to-noise ratio ~10 in 
the core, and ~8 on the arms; its companion, NGC 1510 is detected with only marginal 
significance with a peak signal-to- noise of ~3.5 (Figure [1]). 



Atmospheric emission at 1.1 mm is the dominant noise source in our AzTEC maps. 
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Prev ious images made with the standard AzTEC pipehne (IScott et al. 



2008 



Perera et al. 



20081 ) are optimized for point source recovery and consequently do not preserve extended 
flux in the image. We have developed our own iterative flux recovery technique that 
takes advantage of the fact that the atmospheric noise is not correlated between the 25 
observations and that the atmosphere is transient in time while the astronomical signal 
is stationary. We begi n by using a varia nt of the principal component analysis (PCA) 
technique described in IScott et al.l (|2008[ ) to remove atmosphere from the detector time 
streams. Rather than the aggressive form of flltering described there — which would result 
in a sharp spatial fllter applied to this map — we project just two eigenmodes from the 
detector-detector correlation matrix out of each 15 second span of time-ordered data. This 
relatively weak fllter is comparable in strength to a subtraction of the array average and 
gradient from each detector timestream sample. Once these initial "cleaned" data streams 
are produced we implement the following algorithm: 



1. A coadded image of the fleld is made from the cleaned data streams. 

2. The fluxes from pixels within 24" of pixels with S/N>2.5 are preserved in the cleaned 
image while all other pixels are set to zero. We call this the "current best sky" image 
since it is the best estimate we have of a (noiseless) sky. 

3. The "current best sky" image is recast into the time streams of the detectors and 
subtracted from the original detector time streams, creating a set of residual time 
streams that, in principle, have the same contamination from the atmosphere but less 
true astronomical signal. We call these the "residual observations" . 

4. The "residual observations" are cleaned and mapped. The resulting image is added 
to the "current best sky" image and a new "current best sky" realization is produced 
using the same pixel flux criteria as before. 
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5. This iterative process continues until no new pixels in the realization pass the pixel 
flux criteria. 

6. At this point the severity of atmospheric flltering is increased. That is, another 
detector-detector covariance matrix eigenvector is projected out of the PCA basis for 
the atmospheric cleaning and the iterations continue. 

7. Iterations are complete when the full PCA technique for optimal point source 
detection is realized (in this case cutting 7 eigenm odes in the atmos phere cleaning 
step) and no new pixels pass the pixel flux criteria (IScott et al.l 120081 ). 



One beneflt of this technique is that it results in a set of 25 residual timestreams that 
have been effectively "cleared" of astronomical signal while leaving a good approximation 
of the true atmospheric contamination behind. We take advantage of this by injecting 
various forms of simulated signals into these residual timestreams in order to perform tests 
on the effectiveness of the iterative flux recovery technique. Figure [2] shows the results of 
injecting a flux-scaled versionQ of the MIPS 24 jj,m map of NGC 1512, smoothed with a 
28" Gaussian, into our residual timestreams and then repeating the identical analysis that 
was done for the AzTEC NGC 1512 observations. The residual map is consistent with the 
noise level in the AzTEC NGC 1512 image and shows only minor residual systematics of 
order 1-2 mJy at the locations of the cores of NGC 1510 and NGC 1512. We perform the 
same aperture photometry analysis on the residual map shown in Figure [2] as is done for 
the actual map (see Section 3.3) and flnd that our iterative flux recovery technique misses 

^The flux of the MIPS map has been scaled such that the extended features of the galaxy 
have a flux that matches that of the AzTEC NGC 1512 map. This results in a peak flux in 
the core of the galaxy which is signiflcantly larger than what we actually see at 1.1 mm, so 
this is a robust test of the iterative flux recovery technique. 
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between 0.2% and 12% of the integrated flux depending on the scale-size of the aperture. 
These values are small compared to the total integrated flux and calibration uncertainty 
(17.1% - 128%, see the next section). 



3. DATA ANALYSIS 

3.1. Photometric Apertures 

In addition to measuring the fluxes from the whole galaxies, the infrared and mm maps 
have sufficient angular resolution to allow us to isolate substructures within the larger of 
the two galaxies, NGC 1512. We thus identify three separate subregions in this galaxy: the 
center, and two areas in the spiral arms (Figure [3]). For each of these regions, we define 
circular or ellipsoidal apertures, depending on the structure we are measuring. 

For NGC 1512 as a whole, we define an elliptical aperture (Figure |3]), with major/minor 
axes 491"/287" and centered on the galaxy position defined in NEEo; the aperture's major 
and minor axes are matched to the optical major and minor axes and the position angle 
ex tracted from the KPNO/CTIO i? V7?/ imaging campaign for the SINGS project (Table 1 



of 



Dale et al. 



20071 ). For NGC 1510, which is unresolved in the MIPS and mm images, and 
only marginally resolved in the IRAC bands, we define a single, circular aperture, with 
diameter 108.4", centered again on the NED position of the galaxy. 

The size and location of the photometric apertures for the substructures in NGC 1512 
are more arbitrary in nature, and we define them based on the 1.1 mm map, taking also 



^The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics 
and Space Administration. 
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into account the lower resolution of the 160 /xm image. The central region is defined by a 
circular aperture centered on the NED galaxy center, and with diameter 73.6", delimited 
by the size of the central emission region which is similar in the maps of the two longest 
wavelengths. The location and size of the two apertures along the spiral arms (Figure |3]) are 
chosen to encompass as much as possible of the northeastern and southwestern arm regions, 
while minimizing overlap with the central region. The arm apertures are thus elliptical in 
shape with major/minor axes 174"/78". 



A constant background is determined and subtracted from each Spitzer image, 



using 



Calzetti et al. 



20051). 



the peak value of the background pixels distribution (see description in 
The 1-(T uncertainty of the background band is, in unit of MJy sr~\ [0.008, 0.012, 0.052, 
0.053] for IRAC [3.6, 4.5, 5.8, 8.0] /im, and [0.039, 0.36, 0.48] for MIPS [24, 70, 160] /xm, 
respectively. After background subtraction, photometric values are determined for each of 
the five regions in each band. 



3.2. Photometric Corrections 



Our photometric measurements require a few corrections, due to the presence of: MIPS 
fiux non-linearities, PSF wings outside the apertures, and, for IRAC, scattered light for 
extended sources^ 1.1-mm ASTE/AzTEC data require no aperture corrections. 

MIPS 70 /im pixels with surface brightness above ~66 MJy sr~^ are subject to 
non-linearities that need to be corrected. The central pixel of NGC 1512 is above t 



threshold, and we adopt the non-linearity correctio n formula pub 



(equation 3), which is deduced from data given in ( iGordon et al. 



I is above t 


lis 


(Dale et al. 


2007) 



20071 ). A correction factor 



of 1.011 is derived for the photometry of the whole galaxy, and 1.023 for its central region. 



^http: / / spider.ipac.caltech.edu/staff/jarrett/irac/. 
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The extended source aperture corrections provided by the Spitzer Science Center for 
IRAC gives [0.911, 0.941, 0.798, 0.745] for NGC 1512, [0.928, 0.960, 0.901, 0.791] for its 
central region, [0.919,0.949, 0.866, 0.770] for the north arm, [0.919, 0.949, 0.866, 0.770] 
for the south arm and [0.920, 0.950, 0.869, 0.771] for NGC 1510 at [3.6, 4.5, 5.8, 8.0] /im 
respectively, with an average uncertainty ~10%. 



We adopt 



1512 given in 



the MIPS apertu re corrections ([1.04, 1.05, 1.06] at [24, 70, 160] /im for NGC 



Dale et al. 



(120071 ). For the other apertures we evaluate the corrections by 
measuring the wings of MIPS PSFs for the relevant aperture sizes. The aperture correction 
for the central region of NGC 1512 and for NGC 1510 are [1.06, 1.13, 1.53] and [1.04, 1.08, 
1.34] at [24, 70, 160] /im, respectively. We derive aperture corrections of [1.04, 1.08, 1.28] 
for the arm regions, which are roughly extrapolated from circular apertures inscribed and 
circumscribed to the minor and major axes, respectively, of the aperture ellipses. 



3.3. Uncertainties 



IRAC Calibration un certainties are 5% -10% for 3.6 and 4.5 fim bands and 10%-15% 



for 5.8 and 8.0 /im bands (jPale et al 



120071 ). To be conservative, in this analysis we employ 
the upper limits of those ranges. For MIPS data, 4%, 7% and 12% ar e adopted at 24, 7 



and 160 /im respectively (lEngelbracht et al 



2007 



Gordon et al. 



20071 : 



Stansberry et al 



20071 ). The AzTEC/ASTE 1.1mm image calibration has an uncertainty of ~8^ 



'0. 



Aperture corrections are ~10% uncertain for IRAC channels. For the MIPS bands, the 
aperture correction uncertainties are as small as a few perc ent, and our estim ates of [5%, 



5%, 8%] at [24, 70, 160] /im are close to previous estimates (jPale et al. 



20071 ). 



The background noise levels in the Spitzer images are found by performing a Gaussian 
fit of the pixel-value distributions. For the millimetric image, we test for false detections 
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and fluctuations by performing pliotometry of the background in tlie map using ~8000 
sliglitly displaced consecutive apertures (with the same area as each of the apertures used 
for source photometry) to cover the whole image, but excluding sources and the map edges. 
A Gaussian fit is then applied on the resultant histogram of "background" photometry 
and the dispersion of the background fluctuation is obtained accordingly. This strategy 
results in relative 1-a total uncertainties at 1 mm (errors in absolute flux calibration 
and background estimation added in quadrature) of 17.1% (whole), 20.9% (core), 24.8% 
(northeast arm), 26.7% (southwest arm) for NGC 1512, and 128% for the entire NGC 1510. 
Total uncertainties computed in this fashion are listed in Table |2] for all the IRAC, MIPS 
and AzTEC bands. 



MODEL FITTINGS AND RESULTS 



Draine Sz Lil (120071 ) is used to flt, via 



The silicate-graphite-PAH dust model of 
minimization, the flux densities of the two galaxies and of three sub-galactic regions in all 
eight bands. Parameters in the models are: the PAH mass fraction Qpah, the lower cutoff 
of the starlight intensity distribution, and the fraction of the dust heated by starlight. The 
dust composition is assumed to be a mixture of carbonaceous grains and amorphous silicate 
grains, with size distributions consistent with the observed extinction curve in the local 



Milky Way (Weingartner fc Draine 



20011), along with varying amounts of PAH molecules. 
All charged PAHs (singly or multiply, positively or negatively) of a given size are assumed 
to have the same cross sections, but those of neutral and ionized PAHs are assumed to be 
different. This model employs a flxed shape for the PAH size distribution, a flxed PAH 
ionized fraction, and a flxed spectral shape of the illuminating starlight, which result in 
flxed ratios of PAH emission band strength. The stellar contin uum contributing to the 



mid-infrared emission is approximated by a 5000 K blackbody ( jPraine &: Li 



20071 ) 
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The comparison between data and best-fit models are shown in Figure HI together with 
the resulting model parameters. The SEDs of NGC 1512 and its regions are well represented 
by the models, yielding dust masses Md~2.4xl0'' Mq for the whole galaxy, and in the range 
~2-6xl0^ Mq for the regions within the galaxy (see Table 3), using the distance listed in 
Table 1. NGC 1510 is less well fit by the model, mainly owing to the combination of a high 
70/160 fim ratio (which implies a high effective dust temperature) together with a highly 
uncertain 1.1 mm flux (whose value would imply presence of large amounts of cooler dust, 
see below). Despite the uncertainties, the models still provide a reasonable fit to the NGC 
1510 SED, with a resulting dust mass Md~1.7xl0^ M0, or about 150 times less dust mass 
than its companion galaxy. This difference is reduced to a factor of 45, if we use the upper 
bound of the large uncertainty in NGC 1510's dust mass determination. 

As a test, we have also fit the far-infrared (70 and 160 /im) and mm data to a very 
simple prescription of two modified blackbodies, with temperatures and relative intensities 
derived through minimization (Table 4). The power-law index of the du s t emi ssivity 



is fixed to the value e=2, following the conclusions of both ^Dunne fc EalesI (120011 ) and 



Willmer et al. 



( I2OO9I ) who determine that this value produces better fits to the SEDs of 
the galaxies they consider than lower power-law index values. Overall, this is a very crude 
prescription, with 4 free parameters (warm dust temperature, T„, cool dust temperature, 
Tc, the relative contribution of the two dust components to the observed infrared-mm 
SED, and the total luminosity), thus the fit is not constrained by our three data points 
(70 /im, 160 /im and 1.1 mm). Indeed, we only derive these quantities for comparison 



purposes. However, in order to 
through the physical models of 



guide our p h ysical 



intuition, we use the masses derived 



Draine fc Lil (120071 ) as our "target values" for the simple 



two-temperature prescription. For both galaxies and all regions, we find that the colder 
of the two te mperatures is i n the range Tc~14-16 K (Table 4). As a comparison, within 



the context of 



Draine &: Li 



(I2OO7I ) model, seeing that stronger radiation fields result in 
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hotter dust grains, the lowest dust temperature can be inferred from the lower cutoff of 
the starlight intensity distribution Umin- The best-fit model of the whole NGC 1512 is 
parameterized by f/mm=0.7, corresponding to a lowest dust temperature ~16 K, which 
is consistent with our result of two-temperature fittings as mentioned above. While for 
NGC 1510, Umin=8.0 and therefore a lower temperature limit is found to be 24 K, which 
is somewhat higher than that of the cooler dust component (~15 K) obtained with the 
two-temperature fitting strategy. However, we do not regard this as a discrepancy, given 
that the 1.1 mm data point of NGC 1510 is highly uncertain and that both temperature 
values are still consistent with the fact that active star-forming galaxies appear to have cool 
dust components with temperatures of-order 20 K (see the Discussion section). For the 
warmer dust component, we have a markedly different behavior for the two galaxies: NGC 
1512 and its regions are described by "warm" dust temperatures in the range T„~20-25 K, 
with the higher value being associated with the central region of the galaxy; NGC 1510 
requires a warmer temperature than NGC 1512, T„~36 K, to fit its SED. This is a 
well-known r esult: low-mass ga l axies tend to hav e higher effective dust temperatures than 



large spirals (IHunter et al 



1989 



Dale et al 



20051 ). The total dust masses derived with the 



two -temperature meth od are higher than those derived from the more accurate fits using 



the iDraine fc Lil (l2007l ) models, by a factor 1.6-2.5 in NGC 1512, and by a factor ~9 in 
NGC 1510, or 2.8 if the upper bound of its dust mass is applied (Table 4). Assumptions 
of lower dust emissivity would yield even higher dust m asses from the two-temperature 



method than those derived with the 



Draine &: Lil (120071 ) method (see the Conclusions). 



For NGC 1512 as a whole the warmer dust component includes ~7% of the total dust 
mass in the galaxy, and the dominant contributor by mass is the cooler dust component 
(Table 4). However, in the sub-galactic regions of NGC 1512, the warm dust represents a 
larger fraction of the total dust mass, between ~17% and ~40%. This is consistent with 
the fact that the sub-galactic regions have been selected to encompass actively star forming 
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areas in the galaxy, and are likely biased towards higher warm dust contents than the 
average across NGC 1512. In NGC 1510 the mass-dominant dust component is the cool 
one, with a mass that is about 36 times larger than that of the warmer component; this is 
perhaps not surprising as the higher warm dust temperature for this galaxy relative to NGC 
1512 naturally leads to a much higher emissivity, corresponding to much smaller warm dust 
mass. Finally, we remark again that the results for this galaxy are highly uncertain, due to 
the low-significance of the mm flux determination. 



DISCUSSION 



The availability of mm observations for the nearby, star-forming galaxy pair NGC 
1512/1510 indicates that the observed infrared/mm SED of both galaxies and of the 
three sub-galactic regions in NGC 1512 analyzed in this paper are consistent with the 
presence of a relatively cool dust component (Tc~14-16 K, Figure [3]), if a power- law 
index of the dust emissivity e=2 is assumed. Presence of a non-negligible cool dust 
compone nt with temperature below 20 K ha s already been sugg ested for nearby elliptical 



galaxies ( iHaas et al 



1998 



Klaas et al. 



2001 



Leeuw et al 



2004h . In M 31, a cold dust 



component with a temperature of 16±2 K seems to be present as well ( lOdenwald et al. 



1998 



Genzel &: Cesarsky 



20001 ). Confirming the actual existence of this component in our 



case is difficult, due to the crudeness of the two-temperature models and to the sparseness 
of the SED sampling. In general, data between 160 fim and ~1 mm are not available for 
galaxies, although a few galaxies (including NGC 1512) have been recently observed at 



20091 ). and more 



wavelengths between 250 /im and 500 /im with BLAST (IWiebe et al. 
will become available thanks to the observations with the Herschel Space Telescope. 
Interestingly, even very actively star forming galaxies, like LIRGs and ULIRGs, do seem 



to require some cool dust component, with temperature of-order 20 K, about a factor 
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2-2.5 small er than the temperatu re of the warm dust component that dominates the total 



luminosity (IDunne fc Bales 



20011). 



The observed far-infrared/mm SEDs of the two galaxies cannot be easily fit by a single 
modified blackbody, without requiring that the dust has a power-law index <2: NGC 1512 
requires e~1.3 and NGC 1510 e^l.l, the latter value being at the margin of the acceptable 
range for dust emissivity (l<e<2, Seki & Yamamoto 1980). Less stringent are the fits 
of the sub-galactic regions, where the single-temperature assumption requires e~1.5-1.8, 
marginally consistent with e=2; indeed, these regions have been selected to be actively 
star-forming, so it is not surprising that their infrared emission is dominated by the warm 
dust component. 



Du st emissivity va 



Clouds (lAguirre et al. 



ues 6 <2 have been derived for dust emitting regions in the Magellanic 



20031), assuming a single temperature for the modified blackbody; 



howe yer, two-temperatur e models with e=2 appear to produce equally good fits to the 



data (lAguirre et al. 



20031 ). Furthermore, the multi-wavelength infrared/sub-mm data of 



ULIRGs and SLUGS galaxies seem to indicate that better fits are obtained when multiple- 



temperature modifiec 



(IDunne fc Eales 



2001 



jlackbody emission with emissivity having a power-law index e=2 



Willmer et al. 



20091 ). The pre sence of mu 



Rieke et al 



t iple t empe rature dust 



(120091) and 



Calzetti et al. 



compo nents is supported by the recent conclusions of 
(120091 ). who find evidence for dust self- absorption in the mid-infrared SEDs of LlRGs, 
which then requires re-emission at longer wavelengths. These results provide ground for an 
"onion-peel" scenario, where different dust layers are heated at different temperatures. 

Single-temperature fits yield values around 20-25 K and ~36 K for the warm dust 
components of NGC 1512 (and its sub-regions) and NGC 1510, respectively. This is 
consistent with NGC 1510 being, on average, more actively star-forming than NGC 1512, 
when considering the star formation rate density (Table 1), and also being less metal-rich, 
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which leads to less ability of self-shielding for the dust (jCalzetti et al.ll2000l ). We also notice 
a higher temperature of the warmer dust in the center of NGC 1512 (25 K) compared to 
that of the arm regions of the same galaxy (20-21 K), implying a stronger radiation field in 
the galaxy center where more UV photons are reprocessed by dust grains. 

The combination of a long wavelength baseli ne for the dust em ission (thanks to the 



mm data) and the physically-motivated models (jPraine fc Li 



20071 ) enables us to derive 



accurate dust masses for NGC 1512 and its subregions: typical 1 a uncertainties are in 
the range 19%-25% (Table 3). Less secure is the dust mass of the fainter NGC 1510, 
with a factor ~3 uncertainty, due to its fain t mm emission. We should stress that our 
quoted accuracy is within the context of the iDraine fc Lil ( 120071 ) models. Even with this 



uncertainty, we are in a position of deriving dust-to-gas ratios for the galaxy pair. 

The two interacting galaxies, located in the southern sky, do not have any CO (tracer 
of molecular gas) observations; therefore we will be limite d to the derivati on of Md/Mni 



ratio. Interferometric VLA observations exist for the pair (IThornley 



2006) 



single -dish data from the Parkes Telescope (as part of the HIPASS survey. 



, in addition to 



Koribalski et al. 



20041 ). We use the latter data with low angular resolution to recover the diffuse emission 
in the galaxy pair which is missed by the interferometric data. The Parkes data contain a 
total emission of 259±17 Jy km/s at 21 cm, and by comparing this figure with the total fiux 
(>3cr) contained in the VLA map, we estimate that the interferometric observations miss 
~52% of the total H I gas mass. We assume this "diffuse" component to be homogeneously 
distributed within the area occupied by the two galaxies; this assumption is very simplistic, 
but we expect that its impact on the individual regions (which are dominated by the 
clustered H I detected by the VLA) will not be large, except for the center of NGC 1512 
where the H I flux is weak. We then measure the 21 cm emission within each of the flve 



regions deflned in Table 2 from the VLA image, include the correction for the "diffuse" 
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gas derived above, and derive H I masses using the conventional formula Mhi/M0= 
2.36x10^ l2icmDMpc' where the integrated intensity l2icm is in Jy km/s and the distance in 
Mpc. 

■0 pair, we 



Draine et al. 



When enclosing the entire H I emission from the NGC 1512/NGC 15 
obtain Md/MHi=0.0034±0.0008, consistent with the value 0.0028 found by 
( 120071 ). If only the H I emission inside our photometry aperture for NGC 1512 itself (see 
description in the data analysis section) is taken into account, the Md/Mm ratio becomes 
an order of magnitude higher (0.037). Significantly higher values are obtained for the 
northern (Md/MHi=0.089) and southern (Md/MHi=0.079) arm regions, and for the core 
(Md/MHi=0.13). These values should be, however, considered upper limits to the actual 
dust-to-gas ratios, since we do not have H2 measurements, and molecular gas is likely 
to be an important component in these three actively star-forming regio ns. Based on 



spati a 



2007 



ly-resolved CO s tudies of other massive star-forming galaxies (e.g. 



Leroy et al 



Kennicutt et al. 



20081 ) . we estimate the H2 contribution to the gas content of the 



sub-galactic regions in NGC 1512 to be at least of the same order as the H I itself. Using 
the expected proportionality between t he dust-to-gas ratio and metallicity in galaxies. 



Md/Mgas~0.010(O/H)/(O/H)Mw, from 



Draine et al. 



(l2007l ). we can tentatively compare 



this formula with our derived dust and H I masses. Using the metallicity value in Table 1, 
we find for NGC 1512, Md/Mgas=0.011, which is lower than the Md/Mni value we determine 
the galaxy. We attribute this discrepancy to the absence of molecular gas data. For NG C 
510, we derive Md/MHi=0.0027, and, from its metallicity (Table 1) and 



Draine et al. 



( l2007l )'s formula, an expected Md/Mgas=0.0042. Again, although the two numbers differ by 



about a factor of 2, lack of molecular gas information hampers detailed comparisons. 

While Md/Mni ratios for individual regions and/or the whol e galaxies appear co nsistent 



or larger than what one would expect for Md/Mgas as given by 



Draine et al 



(120071 ). the 
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dust-to-gas ratio of the galaxy pair is lower, by about a factor of 3, than what expected, 
based on the same formulae. This 'deficiency' would be exacerbated by the presence of 
significant molecular gas in the system. We speculate this deficiency to be due to the 
presence of a large fraction of H I not associated with the star-forming disk(s) and/or 
regions within the two galaxies. A solution to this issue will require both higher sensitivity 
H I maps and the availability of CO maps to measure the contribution and distribution of 
the molecular gas. 



6. SUMMARY 

The combination of Spitzer mid/far-infrared images with ground-based ASTE/AzTEC 
1.1 mm data for the galaxy pair NGC 1512/1510 has enabled us to derive accurate 
(~19%-25% uncertainty) dust masses for the large spiral NGC 1512 and three of its 
subregions (the center and the two spiral arms), and constrain the dust mass of the 
low-metallicity dwarf galaxy NGC 1510 to within a factor of 3. The total dust mass of 
NGC 1512 is 2.4x10'^ Mq and it is about 150 times smaller in NGC 1510, in agreement 
with the former being a large spiral galaxy and the latter a small compact dwarf galaxy. 

In both galaxies, when the SEDs are fitted with simplistic two-component modified 
Planck functions, the majority of the dust mass is found to have a relatively cool 
temperature, ~14-16 K, similar to what found for the nearby galaxy M31. Conversely, 
in the central region of NGC 1512 the cool and warm (Tw~25 K) dust contribute close 
to equal mass, reflecting the larger speciflc SFR of the starbursting center relative to the 
galaxy as a whole. The two arm regions display properties that are intermediate between 
the central region and the whole galaxy in terms of ratio between the warm and cool dust 
mass, with values in the range 0.2-0.3. 



The dust-to-gas ratio in NGC 1512 (estimated using only H I, because of lack of 
CO emission data for the pair) is about 3 times below the exp ectation for a galax y with 



Draine et al. 



(120071 ). The 



a metallicity similar to the Milky Way, as already remarked in 
addition of H2 to the gas census would only make the discrepancy worse. Conversely, the 
subregions (center and arms) in the galaxy show high Md/Mni ratios, ~0. 08-0. 09 for the 
arms and ~0.13 in the center, thus about 20 to 30 times larger than the galaxy as a whole, 
and consistent with the lack of H2 information. We speculate that much of the H I included 
in the Md/Mni ratio estimate of the galaxy as a whole is not related to the star-forming 
disk. 



The ASTE project is driven by the Nobeyama Radio Observatory (NRO), a branch of 
National Astronomical Observatory of Japan (NAOJ), in collaboration with the University 
of Chile, and Japanese institutes including University of Tokyo, Nagoya University, Osaka 
Prefecture University, Ibaraki University, and Hokkaido University. A part of this study was 
supported by the MEXT Grant-in- Aid for Specially promoted Research (No. 20001003). KS 
was supported in part through the NASA GSFC Cooperative Agreement NNG04G155A. 
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Table 1: Characteristics of the Galaxies. 





NGC 1512 


NGC 1510 


Morphology^ 


SB(r)ab 


SAO° pec?; H ll BCDG 


fH (km/s) 


898±3 


913±10 


D (Mpc)^ 


10.8±0.8 


11.0±0.8 


12+log(0/H)= 


8.71 


8.31 


V d 
^SFR 


0.001 


0.06 



"Galaxy morphology and heliocentric velocity are from the NASA/IPAC Extragalactic Database (NED) 
''Adopted distances. 
■^Oxygen abundances, from 



Moustakas et al. 



Storchi-Bergmann at al 



(119941) 



(|2009D for NGC 1512 and from 
for NGC 1510. 

"^The star formation rate density, in units of Mq yr^^ kpc^^, from a combination of Ha and 24 /im measure- 



ments for NGC 1512, to account for both the dust obscured and unobscured star formation (Kennicutt et al 



20091 ). and from Ha measurements for NGC 1510 (jCalzetti et al 



20091) . 
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Table 2: Photometric measurements of the galaxy pair NGC 1512/1510. 



A(//m) 


NGC 1512 


NGC 1512c^ 


NGC 1512ne^ 


NGC 1512sw^ 


NGC 1510 


3.6 (mJy) 


349±49 


147±21 


54±8 


52±7 


17±2 


4.5 (mJy) 


220±31 


93±13 


34±5 


33±5 


11±2 


5.7 (mJy) 


244±44 


107±19 


49±9 


45±8 


10±2 


7.9 (mJy) 


394±71 


153±28 


82±15 


79±14 


19±3 


23.7 (mJy) 


447±29 


202±13 


78±5 


73=b5 


136±9 


71.4 (Jy) 


6.24±0.54 


3.24±0.28 


1.03±0.09 


0.81±0.07 


1.04±0.09 


156 (Jy) 


21.3±3.1 


6.65±0.96 


4.99±0.72 


4.68±0.68 


0.98±0.14 


1100 (mJy) 


322±55 


34±7 


62±15 


56±15 


10±12 



tNGC 1512c represents the central region of NGC 1512, while the northeastern/southwestern arms are 
denoted by NGC 1512ne/NGC 1512sw- 
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Table 3: Derived Quantities for NGC 1512 and NGC 1510 



Parameter NGC 1512 NGC 1512c^ NGC 1512ne"^ NGC 1512sw"^ NGC 1510 

Md'' 2.4±0.6xl0^ 2.1±0.4xl06 5.5±1.3xl06 5.7±1.3xl06 1.7±3.6xl0^ 

Mhi,vla^ 4.9x10* 5.4x10^ 4.6x10^ 6.2x10^ 4.7x10^ 

Mmparkes^ 1.7x10* 1.2x10^ 1.5x10^ 1.0x10^ 1.4x10^ 



^NGC 1512c represents the central region of NGC 1512, while the northeastern/southwestern arms are 
denoted by NGC 1512ne/NGC 1512sw. 
Derived dust mass, in units of Mq, from the model fits of [ 



Draine et al 



( 2007 ) 



Extrapolated H i mass, in units of M© (jThornlev 

Mni.Parkcs- 



20061 ). The total H i mass is the sum of Mhi vla and 



-29- 



Table 4: Derived Parameters of Two-Blackbody Fit for NGC 1512 and NGC 1510 



Parameter NGC 1512 NGC 15120^ NGC 1512ne^ NGC 1512sw^ NGC 1510 





24.0 


24.8 


21.4 


20.4 


36.1 


c 


13.8 


16.2 


13.9 


13.8 


15.4 




4.1x10^ 


1.7x10^^ 


1.8x10^ 


2.3x10^ 


4.2x10^ 




5.5x10' 


2.5x10'^ 


8.7 xlO*' 


7.1x10*^ 


1.5x10*^ 



'^NGC 1512c represents the central region of NGC 1512, while the northeastern/southwestern arms are 
denoted by NGC 1512ne/NGC 1512sw 
"Best fit values for the temperatures, and Tc in units of K, from a simple two modified blackbody fit, 

with fixed power-law index of the dust emissivity e=2. 

''Best fit warm and cool dust masses, in units of Mq, associated with and Tc, respectively. 
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Fig. 1.— Mosaic of the images of the galaxy pair NGC 1512/NGC 1510 at 5.8, 8.0, 24, 70, 
160 jJLva and 1.1 mm. North is up. East is left. The field of view of each panel is 8.7'x7.4'. 
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Fig. 2. — (Left) Input MIPS 24 /xm image scaled such that the extended flux matches 
that measured in the AzTEC map. This results in a peak flux at the galaxy center of 
133 mJy. This image is injected into the residual AzTEC timestreams as described in the 
text. (Center) The recovered image using the same FRUITLOOPS algorithm as used to 
produce the AzTEC image of NGC 1512/1510. (Right) The residual image (left image 
minus center image) . Negligible artifacts of 1-2 m Jy from the AzTEC analysis filters remain 
in the cores of the two galaxies due to the bright initial image flux. The faint extended flux 
is recovered to well within the pixel noise of the map. 
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Fig. 3. — The five photometric apertures used in this paper are shown as elhpses/circles, 
superimposed on the 24 /^m map of the galaxy pair. The field of view is 11.4'x8.5'. 
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ngc 1510-whole 

{U) = 13.50, fp^^=0.254 

1pAid=0.004±0.005 
M^=(1.3xlO3±3xl03)(D/Mpc' 



6 

M 

■;^0 




ngc 1512-C( 
j3-2.a, U^,„-3,0, U^^, 
{U)-3.52, f ^,-0.081 
Bt.s=l,0, Ti,,=0,00 
" qp^„= 0,034 ±0.007 
M^={ 1 .Bx 10^±4x^lOa)|D/Mpc; 



^10-'= 




ngc 1512-SW 

(U)=0.62, fp^^=0.031 

qp^f,=0. 039+0. 005 
Mj=(4.9xlO^±l>clO'')(D/Mpc)= 



10 



ngc 1512-corc 

{U)-3.55, f ^,-0.100 

" qp^„=0.034±0.007 
Mj=(l.Bxl0-'±3^^1Q^)|D/Mpc)^^ 

10 , lOf 




ngc 151S-SW 

?=2.0, U™,„=0-6, U^^,= l 
(U)=0.63, f ^^=0.030 
a^g=i.O, Tg.,=0.00 
qp^fi=0-041±0.005 
M^-(4.9xlO''±lxlO'')(D/Mpc)= 



10 



Fig. 4. — The SEDs of NGC 1512, its central region, its northwest and southwest arm 
reg:ions, and its co mpanion NGC 1510. All the SEDs are fitted by the model described in 



Draine fc Lil ( 120071 ). The results of two approaches of the fit are shown in a parallel manner: 
the left panel in each pair allows the parameter /3 to float between 1.5 and 2.5 freely, but it 
is fixed for the right. However, the results sometimes coincide with each other. 



